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SUMMARY 


A  model  of  the  flow  through  a  helicopter  rotor, 
suitable  for  estimating  shed  wake  effects  in  hovering  and 
at  low  tip  speed  ratios,  is  described.  The  method  leads  to 
the  same  result  for  the  aerodynamic  damping  of  blade  bending 
oscillations  as  more  complex  theories.  It  is  shown  that  the 
damping  of  small  natural  oscillations  is  to  a  first  order 
unaffected  by  tip  speed  ratio.  Interharmonic  couplings, 
arising  out  of  the  convection  of  the  shed  vorticity  parallel 
to  the  rotor  disc,  prove  to  be  quite  important,  especially 
if  resonance  occurs.  Also,  in  the  presence  of  strong  inter¬ 
harmonic  coupling,  the  effective  lift  curve  slope  cannot  be 
considered  as  a  universal  parameter,  independent  of  the  har¬ 
monic  content  of  the  blade  motion. 
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I.  INTRODUCTION 


One  of  the  problems  in  the  calculation  of  the 
response  of  rotor  blades  to  the  various  forms  of  excitation 
which  arise,  is  to  assess  the  importance  of  the  'non-steady 
aerodynamic  effects'.  By  this  is  almost  universally  meant 
the  shedding  of  spanwise  vortices  from  the  trailing  edge 
of  a  blade  whose  incidence  is  varying  with  time,  although 
it  really  has  a  wider  connotation.  The  shed  wake  has  an 
effect  on  rotor  blade  oscillations  which  is  an  order  of 
magnitude  greater  than  its  effect  on  the  oscillations  of 
fixed  wings.  Basically  this  is  because  the  wake  does  not 
move  much,  relative  to  the  disc,  in  the  time  taken  for  the 
next  blade  to  occupy  the  position  which  one  blade  has  at 
a  given  instant.  Also  the  blade  natural  frequencies  vary 
with  rotational  speed,  so  that  in  hovering  flight  the  op¬ 
portunity  may  occur  for  shed  vorticity  of  the  appropriate 
sign  and  phase  to  accumulate  underneath,  and  close  to,  the 
rotor,  with  the  result  that  the  aerodynamic  damping  of  flexural 
oscillations  can  become  very  small.  *  But  in  forward 
flight  the  shed  vortices  are  convected  away  from  the  disc, 
and  their  axes  are  no  longer  parallel  to  tin  trailing  edges 
of  successive  blades.  They  therefore  cannot  accumulate  in 
quite  the  same  way  and  their  effect  on  the  rotor  must  be 
expected  to  decrease.  On  the  other  hand  it  is  important  to 
know  just  how  big  is  this  decrease,  for  the  blades  are  then 
in  forced  oscillation  and  should  the  damping  remains  small, 
amplitudes  and  stress  fluctuations  may  be  large. 

But  there  are  very  considerable  difficulties  in 
the  way  of  calculating  non-steady  aerodynamic  effects  in 
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forward  flight.  In  particular  convincing  analytical,  or 
'closed  form'  ,  solutions  are  scarce  because  tractable  models 
of  the  flow  are  hard  to  devise.  The  principal  reasons  are 
the  strong  three-dimensional  character  of  the  flow  -  the 
blade  loading  varies  rapidly  spanwise  and  the  wake  is  at 
least  a  skewed  helix  -  and  the  fact  that  conditions  are  not 
steady  in  either  of  the  two  most  convenient  frames  of  reference. 
A  subsidiary  reason  is  the  sheer  complication  involved,  for 
the  blades  change  their  incidences  at  several  frequencies 
simultaneously  and  there  is  some  degree  of  coupling  between 
frequencies.  Despite  these  difficulties  in  recent  years 
there  have  been  considerable  advances  in  the  technique  of 
calculating  rotor  blade  airloads.  Through  the  combined  use 
of  physical  reasoning,  experiment  ^  9  ^  and  large  computers, 
(5),  (6),  (7)tjiere  kas  been  a  steady  but  striking  improvement 
in  our  understanding  of  the  problems  involved.  But  of  course 
the  various  theories  are  numerical  in  nature  and  involve  large 
and  complex  calculations.  As  such  they  are  not  especially 
suitable  for  establishing  trends  -  except  in  specific  design 
cases.  Moreover,  at  low  tip  speed  ratios,  which  are  of  inter¬ 
est  because  fluctuating  loads  are  large  in  that  region,  there 

are  considerable  problems  in  numerical  analysis  and  some  of 

(8) 

the  basic  approximations  are  invalid.  Piziali  v  ,  for  example, 
has  shown  that  numerical  calculations  of  the  influence  of  the 
shed  wake  are  most  sensitive  to  the  number  of  discrete  vortices 
assumed.  There  is  thus  a  strong  incentive  to  find  some  form 
of  'closed'  solution  so  that  trends  can  be  investigated  and 
a  new  feel  for  the  problem  established. 

This  means  reducing  the  real  system  to  a  model  which 
is  capable  of  precise  and  reasonably  simple  mathematical  form- 
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ulaCion.  It  has  already  been  pointed  out  that  this  is  not 
easy  to  do  and  in  fact  the  conceptual  difficulties  determine 
the  approach,  and  colour  the  procedure,  of  this  paper.  But 
if  the  hovering  case  can  be  regarded  as  solved,  then  in  prin 
ciple  at  any  rate  it  should  be  possible  to  deal  with  flight 
at  very  low  tip-speed  ratios  as  a  short  extension,  or  small 
perturbation,  of  hovering  flight.  Or  conversely  hovering 
may  be  regarded  as  the  limit  of  forward  flight  at  very  low 
tip  speed  ratios.  Thus  we  ought  at  least  to  be  able  to 
estimate  the  rate  at  which  conditions  change  with  tip-speed 
ratio  in  almost  hovering  flight.  For  instance  if  is 

a  measure  of  the  aerodynamic  damping  of  blade  bending  oscil¬ 
lations,  then  it  should  be  possible  to  estimate  the  quantity 
tyi /)/■)/-»<>  •  rhe  purpose  of  this  paper  is,  therefore, 

to  obtain  and  analyse  a  model  of  the  flow  through  a  rotor 
which  provides  a  * closed  form*  solution  for  the  aerodynamic 
damping  the  vicinity  of  JJL  =  0  . 


Now  the  intention  is  to  regard  forward  flight  as 
an  extension  of  the  hovering  case,  so  that  the  first  task 
is  to  establish  a  simple  mathematical  model  for  calculating 
the  effects  of  the  shed  wake  in  hovering.  This  model  must, 
of  course,  be  sufficiently  flexible  for  forward  flight  to 
be  introduced  without  much  additional  analysis.  To  achieve 
this  it  is  necessary  to  make  some  sweeping  assumptions ,  if 
not  in  fact  bare-faced  distortions,  but  these  do  appear  to 
be  justified  by  the  results  -  at  least  for  the  hovering  case 
We  therefore  begin  by  examining  the  important  parameters, 
and  reviewing  previous  closed  solutions  for  the  aerodynamic 
forces  on  rotor  blades  which  are  oscillating  in  hovering 
flight. 
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II.  DEVELOPMENT  OF  THE  FLOW  MODEL 


The  principal  parameters  in  the  problem  of  the 
oscillating  rotor  blade  in  hovering  are  the  number  of  cycles 
of  the  oscillation  which  occur  in  one  revolution  of  the  rotor 
(the  frequency  ratio),  the  pitch  of  the  helical  wake,  the 
number  of  blades  and  the  phase  angle  between  the  oscillations 
of  the  different  blades.  The  frequency  parameter  (reduced 
frequency)  is  of  secondary  importance  although  this  is  only 
because  it  is  small.  Frequency  ratio  is  important  because 
it  determines  how  the  shed  vorticity  is  distributed  over 
the  helical  wake.  A  little  reflection  (or  see  *  ^)), 
will  show  that  if  the  frequency  ratio  is  an  integer  the 
vorticity  will  have  the  same  value  at  corresponding  azimuth 
points  in  every  turn  of  the  helical  wake.  Thus  the  induced 
velocities  of  the  individual  turns  are  cumulative  and  it 
is  this  effect  which  causes  the  aerodynamic  damping  to  be 
reduced.  For  an  oscillation  at  a  given  frequency  ratio  the 
relative  distribution  of  vorticity  in  the  sheets  from  dif¬ 
ferent  blades  depends  upon  the  inter-blade  phase  angle. 

Thus  even  for  oscillations  at  integral  frequency  ratios, 
the  overall  induced  velocities  will  not  be  cumulative  unless 
the  blades  are  oscillating  in  the  appropriate  phase  -  which 
is  a  function  of  the  frequency  ratio  and  the  number  of  blades. 
The  final  magnitude  of  the  induced  velocity  must  also  depend 
upon  how  far  the  wake  vortex  sheets  are  from  the  blades  i.e. , 
upon  the  'pitch'  of  the  wake  which  in  its  turn  depends  upon 
the  inflow  velocity  and  the  number  of  blades. 

The  basic  difficulty  in  the  determination  of  the 
forces  on  a  rotor  blade  is  the  calculation  of  the  induced 
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velocity,  for  the  vortex  sheets  are  spirals  and  their  strength 
is  not  only  initially  unknown  but  varies  both  around  and  across 
the  spiral.  So  far  two  models  of  the  flow  have  been  devised 
which  avoid  or  reduce  this  difficulty.  The  first  is  the  two- 
dimensional  model  of  Jones  ^  and  Loewy  They  argue  that 

since  the  aspect  ratio  of  a  rotor  blade  is  high,  and  most  of 
the  load  must  be  generated  near  the  blade  tips,  then  it  is 
sufficient  to  treat  the  blades  and  their  wakes  as  being  plane 
and  of  infinite  span.  Thus  the  effects  of  the  wake  curvature 
and  spanwise  loading  variations  can  be  ignored  and  each  turn 
of  the  spiral  wake  is  represented  by  a  separate,  plane,  vortex 
sheet  below  the  disc  and  extending  to  infinity  in  front  of  and 
behind  the  blade.  The  sheets  are  spaced  vertically  at  a  dis¬ 
tance  l>/l  which  is  representative  of  the  pitch  of  the  helix, 
and  the  distribution  of  vorticity  in  the  sheets  depends  only 
on  the  frequency  ratio.  One  restriction  on  the  usefulness 
of  this  model  is  that  it  can  be  applied  to  a  multi-blade  rotor, 
only  if  the  inter-blade  phase  angle  is  assumed.  For  our  present 
purpose  the  most  important  prediciton  of  these  theories  is  that, 
at  an  integral  frequency  ratio,  the  two  dimensional  lift  curve 
slope  for  a  vibrating  blade  section  is 


2t 


i. 


k  +  T 


(1) 


i.e. 


the  'lift  deficiency  function'  is  !>/&  ■Mr) 


(5) 

The  second  model  is  due  to  Miller'  '  who  considers 
a  circular  rotor  with  an  infinite  number  of  blades.  In  this 
case  the  frequency  ratio  is  restricted  to  integral  values, 
but  the  influence  of  finite  span  and  wake  curvature  is  included. 
By  considering  a  rotor  which  is  uniformly  loaded  spanwise, 
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so  that  trailing  vorticity  leaves  the  blades  only  at  the  tips, 
the  shed  and  trailing  wakes  are  combined  into  a  circular  cylin¬ 
der  whose  induced  velocities  can  be  calculated.  It  is  assumed 
that  the  variations  in  the  induced  velocity  across  the  chord  of 
a  blade  can  be  neglected  and  lifting  line  theory  is  used  to  calcu¬ 
late  the  aerodynamic  forces.  The  lift  deficiency  function  C 
is  shown  to  be 


so  that  the  number  of  blades  and  the  wake  pitch  appear  through 
the  ratio  •  Here  is  the  solidity  defined  as 

(blade  area/disc  area)  and  )^is  the  ratio  of  the  mean  inflow 
velocity  to  the  tip  speed. 

Equation  (2)  agrees  directly  with  Loewy's  result 
for  C  if  we  put 


The  striking  thing  about  these  two  theories  is 
that  they  lead  to  essentially  the  same  result.  Yet  neither 
theory  is  an  extension  of  the  other;  the  models  used  are  in 
fact  complementary  since  each  contains  features  the  other 
cannot  reproduce.  This  agreement  could  be  a  meaningless 
chance,  for  it  is  possible  that  the  common  prediction  for 
the  lift  deficiency  function  could  be  quite  wrong.  But 
measurements  made  by  Shutler  (reported  in  Jones have 
shown  that  the  predictions  of  Loewy  and  Jones  are  in  good 
agreement  with  experiment,  even  for  the  more  difficult  pro¬ 
blem  of  torsional  oscillations.  Thus  the  evidence  suggests 
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that  the  theoretical  values  for  non-steady  aerodynamic  forces 
are  insensitive  to  the  details  of  the  model  on  which  the  cal¬ 
culations  are  based.  Therefore  why  should  we  not  obtain  a 
very  simple  model  for  the  hovering  case  by  combining  the 
most  (analytically)  convenient  features  of  the  earlier  solu¬ 
tions?  If  this  alternative  representation  yields  the  same 
expression  as  equation  (2) ,  for  the  lift  deficiency  function 
then,  assuming  that  it  has  the  right  flexibility,  it  should 
provide  a  very  good  basis  for  an  extension  to  forward  flight. 

Now  the  simplest  features  of  Miller's  model  are 
the  use  of  lifting  line  theory  and  the  assumption  of  an  infinite 
number  of  blades.  The.  greatest  simplification  introduced  by 
Jones  and  Loewy  is  the  idea  that  the  flow  is  essentially  two 
dimensional.  If  it  is  permissible  to  stretch  these  ideas  to 
the  limit  in  combination  then  we  are  led  to  consider  a  rotor 
which  has  a  very  large  numberof  closely  spaced  blades  each  of 
very  high  aspect  ratio.  For  the  purposes  of  analysis  this 
becomes  a  two-dimensional  actuator  disc,  or  surface  of  dis¬ 
continuity,  exactly  similar  to  that  employed  in  the  most 
elementary  momentum  theories  for  the  propeller.  The  difference 
is  that  since  each  part  of  the  disc  is  to  represent  a  blade 
which  is  oscillating  or  in  some  way  changing  its  incidence 
with  time,  then  we  must  postulate  that  vorticity  is  generated 
at  every  point  of  the  discontinuity. 

A  great  advantage  of  this  model  is  that  the  analysis 
for  it  has  already  been  worked  out  -  or  at  least  the  basic 
procedure  is  established.  Actuator  disc  theory  has  been 
widely  applied  to  the  calculation  of  flow  through  cascades. 

Sears  ^^'has  used  it  to  study  the  rotating  stall  but  the 
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first  real  application  to  non-steady  motion  appears  to  be 
due  to  Whitehead^^  who  used  it  to  calculate  the  forces 
and  moments  on  an  aerofoil  in  an  oscillating  cascade.  Al¬ 
though  we  shall  use  a  different  notation,  and  the  lifting 
line  theory  will  replace  steady  state  cascade  theory,  and 
although  we  find  it  necessary  to  introduce  the  solidity 
to  allow  for  the  finite  number  of  blades  in  a  real  rotor, 
the  analysis  given  below  is  basically  due  to  Whitehead. 

In  fact  much  of  the  stimulus  for  the  present  work  came  from 
Whitehead's  paper. 
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III.  HOVERING  THEORY 


The  actuator  disc  is  shown  in  Figure  1  as  a  dis¬ 
continuity  surface  lying  along  the  axis  of  x,  i.e.,  the 
(infinitesimal)  chord  lines  of  the  airfoils  which  make  up 
the  disc  lie  along  Ox  .  It  is  assumed  that  the  system  of 
co-ordinates  is  rotating  with  the  rotor  blades,  so  that 
these  aerofoils  are  in  a  free  stream  of  velocity  W  where  CJ 
is  the  angular  velocity  of  the  rotor  and  |  is  some  repre¬ 
sentative  radius.  Thus,  relative  to  the  point  at  which  it 
was  shed,  vorticity  moves  parallel  to  Ox  with  velocity  a»l 
and  downwards,  away  from  the  disc,  at  the  mean  inflow  velo¬ 
city  *V=X«o  |  .  Above  the  actuator  disc  there  is  no  vorticity 
and  the  air  moves  uniformly  at  the  mean  inflow  velocity. 

The  flow  model  of  Figure  1  therefore  describes  the  changing 
conditions  on  a  'reference'  blade,  represented  by  the  lift¬ 
ing  line  at  the  origin  0,  as  it  travels  around  the  rotor 
disc.  In  the  analysis  the  actuator  disc  will  be  assumed 
to  extend  to  infinity  in  both  directions,  i.e.  ahead  of 
and  behind  the  reference  blade.  This  is  a  great  analytical 
convenience,  which  has  a  precedent  in  the  work  of  Jones  and 
Loewy,  and  can  be  justified  in  the  following  way.  At  any 
azimuth  position  the  reference  blade  can  be  regarded  as  the 
median  line  -  AB  of  Figure  2  -  of  a  sector  of  the  true  disc. 
On  either  side  of  this  median  the  other  blades,  and  the 
vortices  shed  during  previous  revolutions,  make  an  angle 
with  OB  which,  due  to  the  wake  curvature,  increases  with 
the  azimuth  distance  from  OB.  But  over  some  sector  the 
other  blades  and  the  shed  vortices  can  be  taken  without 
significant  error,  to  lie  parallel  to  the  reference  blade, 
and  the  larger  the  representative  radius  I  ,  the  greater 
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the  arc  on  either  side  of  the  median  -  and  therefore  the 
linear  distance  on  the  model  -  over  which  this  is  sensibly 
correct.  If  this  arc  length  is  many  chords  in  length,  i.e. 
of  the  order  of  the  blade  radius,  then  the  linear  model  may 
be  safely  allowed  to  extend  to  infinity  since  the  remote 
parts,  which  admittedly  are  grossly  unrepresentative,  will 
make  no  significant  contribution  to  conditions  at  the  origin. 


The  equations  governing  the  flow  in  the 
are  those  of  two-dimensional,  inviscid  rotational 
viz,  the  equation  of  continuity 

-^4-^=0 

and  the  dynamical  equations  of  motion 

VF  >-  IS 


r  >x 


x-y  plane 
fluid  motion 

(4) 

(5) 


Vfc 


+  *1 


1  i! 
r  >y 


(6) 


where  Uf  are  the  fluid  velocity  components  parallel 

to  the  x  and  y  axes  respectively,  and  ^  is  the  pressure. 


From  (4)  it  follows  that  there  is  a  stream 
function  X  such  that 

u  -  YL  v  -  _  YY 

u'  ty  >  v» 

The  pressure  p  may  be  eliminated  from  (5) ,  (6) 
to  give  the  vorticity  equation 

)t  +  '  n  *  v'  ^y  ° 


(7) 


(B) 
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where 


(9) 


From  (4),  (7),  (9)  it  follows  that 


y_ x 

lx” 


+ 


(10) 


We  now  assume  that  the  velocities  caused  by  the  shed  vorti- 
city  and  the  blade  motion  are  small  perturbations  of  the 
velocities  4il  and  kl  i.e.  we  put 

a,  =  6ll  +UL  (11) 

V|*'Xwl-«-V  (12) 


so  that, on  neglecting  higher  powers  of  the  perturbation 
velocities,  the  vorticity  equation  becomes 


where 


s 


_  Ivl 


and 


(13) 


For  the  mean,  i.e.  unperturbed,  flow  the  equation  of  con¬ 
tinuity  is  satisfied  by  making  the  same  above 

and  below  the  axis.  But  since  continuity  across  the  actua¬ 
tor  disc  must  also  be  satisfied  by  the  oscillatory  velocity 
components,  it  is  necessary  to  impose  the  additional  condi¬ 
tion  that 


(14) 


In  (14)  the  first  suffix  (zero)  indicates  that 
the  perturbation  velocity  is  evaluated  at  the  disc  (y=o), 
and  in  the  second  suffix,  1  refers  to  the  region  above  the 
disc  and  2  to  the  region  below.  The  distribution  of  vorti¬ 
city  in  the  lower  half  of  the  x-y  plane  -  region  2  -  depends 


11 


upon  the  frequency  ratio  and  the  inter-blade  phase  angle. 

In  the  two-dimensional  model  with  discrete  wakes  these 
parameters  can  be  regarded  as  separate,  but  in  the  present 
model  they  cannot  be  distinguished  because  the  number  of 
blades  is  supposedly  infinite.  Moreover  the  concept  of  a 
phase  angle  is  not  so  appropriate  to  the  linear  actuator 
disc  as  it  is  to  the  real  circular  rotor.  A  more  suitable 
idea  is  to  assume  that  a  wave  of  disturbance  is  propagating 
along  the  surface  of  discontinuity.  For  suppose  the  vortic  ity 

shed  by  the  reference  blade  to  be  of  the 

form 

S(o,0,-t)  =  S^c.*  f+  +  Sfs  *i*i  j>t  (15) 

Then  the  vorticity  shed  by  a  blade  spaced  an  angle  oC 
from  the  reference  blade  will  be 

"Ss:  (16> 

where  p£  is  a  phase  angle  -  or  more  precisely  f  is  a 
time  delay  -  which  in  general  will  be  a  function  of  cL  and 
p  •  But  on  the  linear  actuator  disc  •^-becomes  a  distance, 

Xo  say,  so  that  the  time  delay  must  be  transformed 

into  a  dependence  upon  X0  .  This  is  easily  done  if  we 
put  £  *»  y./u,  where  I  is  a  characteristic  speed 

which  depends  upon  the  frequency  ratio  and  the  inter-blade 
phase  angle.  Then 

&J+ V,B  (17) 

which  repeats  itself  at  intervals  of  time  of  and  inter¬ 

vals  of  of  lirTZ/f  i.e.  it  is  the  equation  of  a  wave 
travelling  along  Ox.  with  speed  . 
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The  vorticity  at  an  arbitrary  point  R(x,y) 

can  be  related  to  the  way  in  which  vorticity  is  shed  from 
the  disc,  since  (see  Figure  l)>the  vorticity  at  ^ left 
at  an  earlier  time 


Thus 


v  )Wk ; 


But 

Hence 


*(*lY f/0»  C#S|> 

+  5fS  «*p(t-  v<r-  'fe) 
x.  *  x-uAy/v 

sfo,y/t)=  T*r*~  v"  0“  wj) 

4  5f*  si*  "tr  =  W 


(18) 

(19) 


(20) 


Equation  (20)  show?  that  the  vorticity  distribution  in 
the  direction  depends  upon  the  value  of  4>-4^i.e.  upon  the 
absolute  propagation  speed  of  the  wave  which  determines  the 
inter-blade  phase  angle.  We  can,  of  course,  proceed  with 
the  analysis  without  specifying  the  ratio  ^  but 
there  is  a  considerable  simplification  if  at  this  stage  we 
introduce  the  value  appropriate  to  the  blade  oscillations 
which  occur  in  forward  flight.  In  forward  flight,  conditions 
at  each  azimuth  are  determined  only  the  cyclic  variations 
in  relative  wind  speed;  they  are  certainly  independent  of 
whichever  blade  happens  to  be  passing  through  that  azimuth. 

Thus  the  various  harmonics  of  the  blade  motion  must  adjust 

their  phase  so  that  at  each  azimuth  conditions  on  any  one  blade 

are  the  same  as  those  which  obtain  on  any  other  blade  at  that  azimuth. 

To  see  how  this  is  achieved  consider  any  two  blades  spaced  an 

angle aCapart  in  azimuth,  and  let  some  quality  <j^(life ,  blade 

motion,  etc.)  associated  with  the  leading  blade  be  of  the 

form  Qj 3  Jdt  •  Then  the  same  quantity  for  the  second 

blade  must  be  of  the  fojrm  where  A  is  the 
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phase  angle  which  we  must  determine.  But  when  4 the 
second  blade  will  replace  the  first  at  any  chosen  azimuth,  and 
the  variable  (J*  is  now  equal  to  A**  ft/*)  .  Thus 

and  .  The  angle  cC  is  essentially  a 

distance-measured  around  the  azimuth,  so  that  the  distribution  of 
inter-blade  phase  angle  is  equivalent  to  a  wave  travelling  around 
the  disc  at  the  rotational  speed.  Alternatively,  we  note  that 
for  the  sort  of  oscillations  which  occur  in  forward  flight, 
an  observer  outside  the  rotor  disc,  but  moving  with  the 
velocity  of  the  hub,  sees  that  any  variable  such  as  G(  varies 
from  azimuth  to  azimuth,  but  that  its  value  at  a  given  azimuth 
is  independent  of  time.  To  this  observer  the  distribution  of 
around  the  disc  takes  the  form  of  a  stationary  wave,  and 
oscillations  of  this  type  can  be  represented  on  the  actuator 
disc  model  by  putting  4J»AJt  .  This  substitution  considerably 
simplifies  the  analysis  since  it  makes  the  vorticity  distri¬ 
bution  independent  of  y  .  At  the  same  time  it  illustrates 
the  connection  between  interblade  phase  angle  and  frequency 
ratio,  since  the  same  independence  of  y  is  observed  if  the 
blades  of  a  hovering  rotor  are  oscillating,  in  the  appropri¬ 
ate  phase  relation,  at  an  integral  number  of  cycles  per  re¬ 
volution.  The  method  can,  of  course,  accommodate  arbitrary 
combinations  of  frequency  ratios  and  phase  angles  but  the 
proper  value  of  must  then  be  worked  out  separately  for 
each  case. 


The  usual  application  of  actuator  disc  theory  is 
not  to  the  calculation  of  the  forces  on  individual  blades, 
but  to  the  determination  of  the  forces  exerted  by  the  rotor 
as  a  whole.  (The  reactions  of  the  fluid  which  balance  these 
forces  can  be  calculated  very  easily  by  applying  the  momen¬ 
tum  and  energy  theorems  of  fluid  mechanics.)  Miller  ^ 
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has  in  fact  used  this  approach  to  show  that  the  lift  defi¬ 
ciency  function  for  the  response  of  a  complete  rotor  to 
changes  in  collective  pitch  is  given  by  (2),  This  must  be 
one  of  the  first  applications  of  simple  momentum  theory  to 
non-steady  motion,  (although  its  use  for  quasi-steady  pro¬ 
blems  has  been  justified  (12)),  and  it  provides  an  impor¬ 
tant  extension  to  classical  propeller  theory. 

But  when  the  incidence  changes  are  periodic  this 
approach  is  not  entirely  adequate,  since  only  those  fre¬ 
quencies  which  are  simultaneously  integral  multiples  of 
the  rotational  speed,  and  of  the  number  of  blades,  can 
contribute  to  the  resultant  force  on  ♦'he  rotor.  It  is  then 
necessary  to  use  a  theory  which  takes  into  account  the  con¬ 
ditions  experienced  by  individual  blades  as  they  rotate.  The 
simplest  of  this  type,  and  the  one  we  shall  adopt  here,  is 
the  lifting  line  theory.  The  validity  of  this  theory  is 
questionable,  particularly  for  large  frequency  parameters 
or  if  the  'in-phase'  component  of  the  aerodynamic  force  is 
required,  but  these  restrictions  are  quite  compatible  with 
the  other  assumptions  which  have  been  made. 

Let  the  circulation  about  a  blade  at  the  repre¬ 
sentative  radius  be  f(;  then  the  lift  per  unit  span  L  at 
that  radius  is 


L  =y»“>l 

(21) 

But 

K  •=  k  CL<ol 

(22) 

Where 

i+w  i 
cot  J 

(23) 

Here  «^is  the  mean  incidence,  %  is  the  upward  displacement 
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of  the  blade  section  and  W  is  the  induced  velocity  -  mea¬ 
sured  positive  downwards-of  the  shed  vortex  system  -  all 


evaluated  at  the  representative  radius*  Considering  only 
the  time  dependent  part  of  K(*Kt) 


(24) 


Note  that  (24)  does  not  contain  the  relative  wind  speed,  so 
that  the  same  relation  between  l^and  i  holds  in  forward 
flight.  Now  in  an  aerodynamic  calculation  must  be  re¬ 
garded  as  known  since  it  is  just  as  much  a  function  of  the 
structural,  as  of  the  aerodynamic,  properties  of  the  blade. 
Thus  to  calculate  the  lift  as  a  function  of  t  it  is  necessary 
to  determine  the  relation  between  and  W  .  But  the 
intention  is  to  use  actuator  disc  theory  to  evaluate  W 
i.e.,  W  is  calculated  on  the  assumption  that  the  number 
of  blades  is  infinite,  whereas  l^is  the  circulation  around 
a  single  blade  with  a  finite  level  of  loading.  We  have 
therefore  to  assume  a  relationship  between  the  circulation 
about  a  finite  blade  and  the  distribution  of  vorticity  in 
the  surface  of  actuator  disc. 


Miller  ^overcomes  this  difficulty  by  in  effect 
assuming  the  blade  loading  to  be  distributed  ever  part  of 
the  circular  disc  and  the  shed  vorticity  to  be  distributed 
over  the  corresponding  part  of  the  wake.  This  approach  is 
not,  however,  suitable  to  the  linear  model  in  which  the  fre¬ 
quencies  and  phase  angles  are  arbitrary.  Instead  we  assume 
the  loading  at  a  point  on  the  actuator  disc  to  be  the  same 
as  the  mean  density  of  loading  on  an  actual 

blade  at  the  corresponding  azimuth.  This  of  course  gives  values 
for  the  induced  velocity  w  which  are  far  too  high,  because 
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the  rotor  is  being  assumed  to  be  much  too  'solid'.  But  since 
the  system  is  linear  the  increase  of  scale  is  the  only  dis¬ 
tortion  introduced  i.e.,  no  phase  shifts  or  additional  fre¬ 
quencies  are  created  by  assuming  every  point  of  the  actua¬ 
tor  disc  to  be  much  more  heavily  loaded  than  it  really  is. 
Therefore,  when  we  come  to  calculate  the  induced  velocity 
at  an  actual  blade,  it  is  merely  necessary  to  multiply  the 
actuator  disc  value  for  this  velocity  -  say  -  by  some 
factor  which  allows  for  the  finite  number  of  blades  in  the 
real  rotor.  The  obvious  factor  to  use  is  the  solidity 
i.e.  W  *  T  Wl  (25) 

where  XT  is  the  solidity  appropriate  to  the  linear  actuator 
disc  theory,  viz.  y-  which  is  the  fraction  of  the  represen¬ 
tative  circumference  occupied  by  blade  cross-sections. 

To  calculate  the  loading  on  the  reference  blade 
we  have  to  determine  the  pressure  difference  across  it,  and 
then  relate  this  pressure  difference  to  the  circulation 
about  the  blade  and  hence  to  the  vorticity  which  is  shed 
into  the  wake. 

From  the  linearised  dynamical  equations 

f> fyok>U  =  -y>  f (it  u  \  iy 

so  that  the  pressure  difference  aP  between  the  lower  and 
the  upper  surfaces  at  a  point  on  the  reference  blade  is 

aP=. 

Now  the  circulation  around  the  reference  blade  is 

K*. = 

so  that  in  accordance  with  the  lifting  line  theory  the 
force  per  unit  span  on  it  is 


(26) 


(27) 
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(28) 


yOO)  l  J  ^ 


But 


L 


is  also  equal  to  the  integral  of  the  pressure  dif¬ 


ference  across  the  chord 


i.e, 


wl  f  far  0<^  -/>& 

cUrJ.  '  tW 


(29) 


Reversing  the  order  of  integration  in  the  last  term  of  (29) 

{jytolft**  O  J  x y  JJx  (&  Ato)  Odx  ( 30  > 

On  comparing  (28)  and  (30)  it  follows  that  we  must  have 


But  this  is  precisely  true,  since  (31)  merely 
expresses  the  fact  that  the  changes  in  the  circulation  about 
a  blade  are  convected  away,  with  velocity  V*  normal  to  the 
disc,  as  vortices  of  opposite  sign,  and  this  idea  has  al¬ 
ready  been  used  to  establish  the  distribution  of  vorticity 
in  the  wake  -  see  equations  (19),  (20). 


An  alternative  form  of  (31)  is 

“'AJ  5.  <32> 

where  5^  is  the  vorticity  shed  at  the  origin  i.e.,  at  the 
reference  blade. 

Whitehead^^  arrives  at  the  same  result  by  im¬ 
posing  the  condition  that  the  difference  in  total  head 
across  the  actuator  disc  is  independent  of  X  •  Since 
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Pyulu.  is  the  linearised  total  head  this  assumption  can  be 


stated  as 


^•(Pi  +^»U)*o 


(33) 


and  (32)  follows  at  once  from  (26) • 


From  (21)  and  (24)  lifting  line  theory  gives  for 
the  fluctuating  part  of  the  lift 


\_  =  _  #  &  (i+Vi) 


(34) 


where  and  W*  is  the  normal  component  of  the  in¬ 

duced  velocity  at  the  reference  blade  i.e. 


from  (14) 


(33) 


Substituting  (35)  in  (34) ,  and  equating  the  two 
expressions  for  the  lift  we  have 


/  -  -aUi-vT  v„ ) 


(36) 


'«W 


We  next  assume  that  the  integral  across  the  chord 
is  equal  to  the  value  of  at  the  reference  blade 

multiplied  by  the  chord,  i.e.,  it  is  assumed  that  the  pressure 
difference  is  constant  across  a  chord  length.  This  seems  to 
be  a  legitimate  assumption  provided  that  the  frequency  para¬ 
meter  is  small. 


Hence 


So  far  we  have  three  equations,  (14),  (32),  (37) 


(37) 
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to  be  solved,  in  terms  of  the  velocity  £  ,  for  the  four 

unknowns  .  v/.  V  A  fourth  relation  can  be  ob- 

J  J  *0* 

tained  from  (10) • 


To  solve  this  partial  differential  equation  the 
variables  are  first  put  into  non-dimensional  form  by  means 
of  the  transformations 


(38) 


Then  substituting  equation  (20)  for  the  vorticity  in  (10) 
we  obtain 

The  solution  for  p”  consists  of  a  complementary  function 
plus  a  particular  integral.  Above  the  disc,  where  there  is 
no  vorticity,  the  particular  integral  is  zero.  Below  the 
disc  the  vorticity  causes  the  particular  integral  of  (39) 
to  vary  periodically  with  ^  and  the  complementary  functions 
above  and  below  the  disc  must  posess  this  same  periodicity. 
This  if  is  the  non-dimensional  stream  function  in  the 
upper  half  of  the  plane 


3 


(39) 


(AO) 


where  the  are  arbitrary  constants  which  are  to 

be  determined.  Either  the  positive  or  negative  exponent 
is  possible  but  the  only  permissible  solution  is  one  which 
remains  finite  as  tends  to  infinity.  Therefore,  since 

is  positive  above  the  disc  we  take 
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\ 


(41) 


The  complementary  function  for  the  lower  half 
of  the  plane  has  a  similar  form,  but  the  positive 

exponent  must  be  used  in  that  region.  Thus  the  complete 
non-dimensional  stream  function  below  the  disc  is 

n 6f-  f  ^  )J  e^Lrc^l 

4-  &«■  t«n  (l-  f)  +  i’lK-A 


(42) 


where  '  c^.s*  are  arbitrary  as  before.  (42)  is  easily 
verified  by  direct  substitution  in  (39). 

Now  the  normal  components  of  the  perturbation 
velocity  must  be  continuous  across  the  actuator  disc,  i.e., 
which  becomes 


(S.  •  (i). 


(43) 


where  the  notation 


means  that  the  derivative  is 


evaluated  along  the  x-axis. 


Making  use  of  the  expressions  (41) ,  (42)  for 
and  and  equating  the  coefficients  of 

Ml  separately  to  zero,  equation  (43)  becomes 

Cw*  Ctw*  <“> 

C  _  C  X-  J**#  tn 


(45) 
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la  non-dimensional  form,  Che  relation  (32)  between 
the  bound  and  shed  vortices  becomes 


which,  on  substituting  for  ^  ^  and  and  equati 
efficients  of  cofrVy^K”^  ’ - 


ing  co¬ 


reduces  to 


^  m_=  ^  ^w«  / 

S(V|  *  \  i  /'A 


Final!}’  the  non-dimensional  form  of  (37)  ,  which 
is  a  form  of  the  relation  between  the  bound  vorticity  and 
the  circulation  around  the  blades  becomes 

GHi-Hf-©} 

To  solve  (49)  we  have  to  assume  that  j  is  periodic  with 
frequency  ratio  X. 

Hence  we  put 

=t  =1  *m+t} 

Substituting  for  ^  and  in  (49)  ,  and  equating  the  co¬ 
efficients  of  Ctfn^^wJj^we  get 

C„+  C„  '  If*,.  +T 


(46) 


(47) 

(48) 


(49) 


(30) 


(51) 

(52) 


22 


Thus, in  terms  of  the  blade  vibrational  velocity,  we  have  six 
equations,  (44) ,  (45),  (47),  (48),  (51),  (52)  for  the  six 
unknowns  Eliminating  the  arbitrary 


constants  we  obtain 


jr]  ^  ’>*-*  0 

The  fluctuating  lift  per  unit  span  on  the  reference 


blade  is 

L  =yo“l«-6ul-*U. bit  ™ 

cash|^  (56) 


Solving  (53),  (54)  for  ^^{j^and  substituting  in  (56) 
have  finally  \ l  , 

=~TUyiV^n 


we 


where 


Now  the  expression 


ai-uii  .  h 


(57) 


(58) 


is  the  non-dimensional  form  of  the  blade  vibrational  velocity, 
so  that  the  lift  force  is  out-of-phase  with  this  velocity 
by  an  amount  which  depends  upon  .  But  usually 

very  small,  and  if  we  neglect  terms  in  £  and  ignore  ^ 
a  first  approximation  to  (57)  is 
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Thus  the  lift  primarily  is  proportional,  but  of  opposite 
sign,  to  the  vibrational  velocity  i.e.,  it  exerts  a  damping 
force  on  the  oscillation.  This  approximate  result  can  also 
be  derived  more  quickly  by  eliminating  the  unimportant  cross¬ 
couplings  from  (53) ,  and  from  (54) . 

The  more  complete  formula  (57)  shows  that  there 
is  in  fact  a  component  of  force  in  phase  with  the  displace¬ 
ment  i.e.,  the  aerodynamic  forces  also  supply  a  stiffness. 
However,  this  term  can  be  ignored  without  difficulty  since 
it  will  make  a  negligible  change  to  the  total  stiffness  of 
r.  blade,  and  the  lifting  line  theory  is  in  any  case  not 
accurate  for  the  'in-phase'  components. 

(a/r) 

The  term  is  the  lift  deficiency  func¬ 
tion  for  the  oscillations.  Observing  that  {T  this  re- 

A  *  (l) 

duces  to  the  same  formula  for  ^  as  was  obtained  by  Loewyv  ' 

and  Millcr^^. 

We  have  thus  developed  a  method  of  calculating 
the  lift  deficiency  function  which  gives  the  same  results  as, 
and  is  easier  in  use  than,  earlier  theories.  The  next  step 
is  to  extend  the  calculation  to  forward  flight  and  this  means 
stretching  the  model  to  its  utmost,  for  we  are  restricted  to 
dealing  with  time-varying  problems  in  two-dimensions.  Some 
indulgence  has  therefore  to  be  granted  when  reading  the 
following  sections,  but  it  is  hoped  that  the  approach  will 


24 


be  taken  seriously  for  the  main  purpose  is  not  so  much  to 
provide  numerical  results,  as  to  establish  a  new  vehicle, 
or  framework,  of  understanding  for  the  forward  flight  prob¬ 
lem.  Despite  its  limitations  the  theory  bears  many  of  the 
hallmarks  of  a  forward  flight  analysis  -  complications,  the 
existence  of  many  frequencies,  differential  equations  with 
periodic  coefficients,  inter-harmonic  couplings,  etc.,  some 
of  which  have  not  been  considered  before.  Furthermore  a 
perturbation  procedure  is  used  -  not  for  the  first  time  in 
rotor  analysis  -  but  in  an  aerodynamic,  as  distinct  from  a 
dynamic,  calculation.  The  use  of  these  techniques  as  a  means 
of  simplifying  rotor  analysis  has  recently  been  advocated 

(13) 

by  Young v  '  and  it  is  hoped  that  this  example  of  their  use 
will  provide  a  basis  for  future  development.  For  example  it 
may  be  that  the  best  approach  to  the  calculation  of  rotor 
loadings  and  noise  at  low  tip  speed  ratios  is  to  expand  the 
wake  geometry  in  powers  of  the  tip-speed  ratio  -  as  well  as 
other  parameters. 
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IV.  EXTENSION  TO  FORWARD  FLIGHT 
4.  I  Forward  Flight  Model 

The  principal  differences  between  hovering  and 
forward  flight  are  that  the  inflow  velocity  and  the  tan¬ 
gential  component  of  the  relative  wind  speed  vary  period¬ 
ically  as  a  blade  rotates,  and  that  the  blades  are  in  forced 
oscillation.  The  relative  wind  speed  variation  has  a  period 
of  once  per  revolution,  but  the  induced  velocity  and  the 
blade  motion  may  have  significant  content  at  all  frequencies 
up  to  ten  per  revolution. 


The  induced  velocity  variations  are  easily  intro¬ 
duced  into  the  model  of  Section  2,  because  the  lifting  line 
theory  does  not  recognise  a  difference  between  the  velocity 
of  the  air  and  the  velocity  of  a  blade  normal  to  the  disc. 
Hence  we  merely  replace  the  velocities  jtx  of  the 

hovering  analysis  by  the  relative  air  velocities  .SJX.W 
where  AJW  is  the  (non-dimensional)^*1  harmonic 
component  of  the  induced  velocity.  The  summation  over  n is 
necessary  because  so  many  harmonics  are  present  and,  as  we 
shall  see  below,  there  is  some  mutual  interference  between 
them.  For  the  purposes  of  the  present  analysis  we  must  as¬ 
sume  that  the  are  known  from  some  other  source,  and 

that  they  are  not  affected  by  any  blade  motion  which  they 
produce.  This  is  probably  a  legitimate  assumption  since, as 
Miller  has  pointed  out^^  ,  these  induced  velocity  varia¬ 
tions  are  primarily  due  to  the  shape  and  mean  strength  of 
the  trailing  vortex  system,  parameters  which  are  fixed  by 
performance  considerations. 


But  though  the  induced  velocity  is  assumed  to 
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cause  a  variation  ot  incidence  around  the  disc,  we  shall  as¬ 
sume  in  the  theory  that  shed  vorticity  is  convected  away 
from  all  points  of  the  disc  at  the  mean  induced  velocity. 

This  assumption  of  a  ’rigid'  wake  is  certainly  not  valid,  but 
not  much  is  yet  known  about  the  movement  of  vortices  after 
they  have  left  the  blades  (for  a  review  of  some  of  the  data 
see  (15)).  When  this  behaviour  is  better  understood  it  will 
be  a  straight-forward,  if  complicated,  matter  to  include  iL 
in  the  actuator  disc  theory.  Until  that  time,  however,  the 
rigia  wake  approximation  is  probably  one  of  the  weakest  parts 
of  the  present  analysis. 

The  periodic  wind  speed  variation  has  two  effects. 
The  first,  which  through  the  lifting  line  theory  is  easily 
absorbed  into  the  actuator  disc  analysis,  is  that  there  are 
periodic  lift  changes  even  if  the  circulation  around  a  blade 
is  constant  with  time.  The  second  effect  is  that  because 
the  wind  direction  is  changing,  the  associated  shed  vorti¬ 
city  is  convected  away  from  the  blade  in  a  direction,  and 
at  a  speed,  which  depends  upon  the  position  of  the  blade  in 
azimuth.  Thus  if  a  blade  is  on  the  advancing  side  shed  vor¬ 
ticity  is  convected  away,  at  right  angles  to  the  blade,  at 
a  speed  which  is  greater  than  that  due  to  rotation  alone, 
and  if  the  blade  is  on  the  retreating  side  the  convection 
is  reduced.  For  a  blade  lying  fore  and  aft  the  convection 
velocity  normal  to  the  blade  is  unaffected,  but  shed  vorti¬ 
city  will  move  parallel  to  the  blade  towards  the  rear,  and 
perhaps  out  of,  the  disc.  Accounting  for  this  convection 
provides  the  main  conceptual  difficulty  of  our  model. 

In  forward  flight,  therefore,  conditions  are  not 
steady  even  in  a  rotating  frame  of  reference,  and  at  a  given 
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instant  these  conditions  vary  from  blade  to  blade.  It  follows 
that  to  represent  forward  flight  on  the  actuator  disc  model, 
either  these  features  must  be  incorporated  in  some  way  or 
they  must  be  shown  to  be  unimportant.  Now  obviously  the 
variation  of  relative  wind  speed  with  azimuth  cannot  be 
ignored,  since  this  is  one  of  the  essential  differences 
between  forward  and  hovering  flight,  but  it  is  probably  quite 
reasonable  to  neglect  its  change  from  blade  to  blade  of  the 
actuator  disc  model.  In  fact  to  be  consistent  with  the  as¬ 
sumptions  of  the  hovering  theory,  this  change  should  be  ne¬ 
glected.  The  reason  is  that  the  actuator  disc  model  sup¬ 
posedly  represents  that  part  of  the  flow  which  is  within 
a  few  chords  length  of  the  reference  blade.  (It  will  be 
remembered  that  -  see  Section  2  above  -  the  linear  model  is 
allowed  to  extend  to  infinity  only  as  an  analytical  conven¬ 
ience.)  At  the  representative  radius  a  few  chords  length 
corresponds  to  a  small  arc  of  azimuth,  so  that  within  the 
compass  of  the  model  the  variation  in  relative  velocity  is 
small  and  can  be  neglected,  i.e.  it  is  assumed  that  the 
relative  wind  speed  at  every  point  of  the  actuator  disc  is 
the  same  as  it  is  on  the  reference  blade,  but  of  course  this 
relative  speed  varies  with  time.  Exactly  the  same  sort  of 
argument  was  used  by  Jones and  Loewy^'to  justify  the 
neglect  of  the  curvature  of  the  wake  in  the  hovering  theory. 
These  approximations  would  not  be  valid  if  the  tangential 
velocity  varied  many  times  per  revolution.  Nor,  by  analogy, 
is  the  argument  justifying  the  neglect  of  wake  curvature  in 
hovering  valid  if  the  frequency  of  incidence  change  is  many 
cycles  per  revolution,  for  then  the  zone  of  influence  of 
the  shed  vorticity  cannot  be  considered  as  limited  to  a  few 
chords  length  on  either  side  of  the  reference  blade.  But 
it  is  clear  from  the  sketch  of  Figure  2  that  the  reduction 
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in  induced  effects  due  to  wake  curvature,  and  again  by  analogy, 
the  variation  in  tangential  wind  speed  over  quite  a  large  por¬ 
tion  of  the  disc,  depend  only  on  the  semi-included  angle  sub¬ 
tended  at  the  hub  by  the  arc  at  the  representative  radius. 

These  reductions  and  variations  will  be  small  for  a  semi¬ 
angle  of  up  to  30°  i.e.  a  total  arc  length  of  about  one  blade 
radius  which  is  sufficiently  large  to  ensure  that  all  the 
important  induced  effects  are  adequately  represented  by  a 
linear  model. 

The  general  convection  of  vorticity  behind,  and 
tc  the  rear  of,  the  disc  means  that  a  blade  may  never  again 
pass  close  to  a  vortex  which  it  shed  in  an  earlier  revolu¬ 
tion,  or  that  if  it  does  so  the  shed  vortex  will  probably 
be  no  longer  parallel  to  the  blade.  This  point  is  illustrated 
in  Figure  2.  AA  represents  a  vortex  shed  from  a  blade  in  the 
position  OC.  . ssuming  the  rotor  to  be  fixed  in  space  and 
air  to  be  flowing  past  it  at  the  forward  speed  TJj.  the  vor¬ 
tex  AA  will  move  to  the  position  A' A'  in  the  time  it  takes 
OC  -  or  some  other  blade  -  to  reach  the  position  OB.  But 
A'A'  is  substantially  at  right  angles  to  OB  so  that  its 
influence  on  OB  will  be  very  different  from  that  of  AA  on  OC. 

Thus  in  forward  flight  the  induced  effects  of  a 
vortex  vary  with  time  because  the  angle  which  it  makes  quite 
different  angles  with  successive  blades.  So  far  no  method 
has  been  devised  of  including  this  feature  in  the  two-di¬ 
mensional  actuator  disc  model,  but  it  is  thought  that  the 
error  in  omitting  completely  will  not  be  too  great.  The 
reason  is  that  the  induced  effect  of  a  particular  vortex 
diminishes  with  time  anyway  since  the  downwash  moves  it 
away  from  the  disc.  In  reality,  particularly  at  low  tip 
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speed  ratios,  the  vortex  A' A'  would  be  well  below  the  disc 
at  the  instant  shown  and  its  absolute  contribution  to  the 
induced  velocity  would  be  small.  This  particular  considera¬ 
tion  does,  however,  suggest  that  this  model  of  the  flow  is 
only  likely  to  be  applicable  to  tip  speed  ratios  -JL- 
which  are  less  than,  or  about  equal  to,  the  mean  inflow 
coefficient  . 


In  practice  the  most  likely  source  of  error  will 
be  in  the  estimation  of  the  effects  of  vorticity  shed  when 
the  blades  lie  more  or  less  along  the  line  of  flight.  Since 
the  number  of  blades  is  limited,  a  finite  time  must  elapse 
before  the  next  blade  enters  the  field  of  influence  of  such 
a  shed  vortex.  During  that  time  the  vortex  may  have  moved, 
parallel  to  itself,  perhaps  out  of  the  disc  or  out  of  the 
region  of  substantially  two-dimensional  flow,  and  it  may 
never  affect  the  next  blade  at  all.  (Of  course,  if  the  blade 
were  really  of  infinite  radius  and  the  flow  strictly  two-di¬ 
mensional,  moving  the  vortex  parallel  to  itself  would  make 
no  difference.)  It  is  possible  then  that  the  theory  will 
over-estimate  the  influence  of  the  shed  wake  in  these  regions, 
but  three  dimensional  and  non-rigid  wake  effects  are  likely 
to  be  of  at  least  equal  importance,  so  there  is  little  point 
in  attempting  to  refine  the  theory  at  this  stage. 

Since  this  is  a  two-dimensional  model,  spanwise 
loading  variations  cannot,  be  included  and  the  effects  of 
wake  curvature,  distortion  and  inclination  due  to  forward 
speed  have  already  been  rejected  or  ignored.  The  final 
actuator  disc  model  for  forward  flight  is  therefore  the 
same  as  that  for  the  hovering  case  except  that  the  relative 
tangential  and  induced  velocity  components  are  allowed  to 
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vary  with  time.  Thus  only  these  'three-dimensional'  effects 
which  can  be  expressed  as  time  variations  are  included, 
and  the  best  justification  for  this  whole  elementary  approach 
is  that  it  is  really  only  meant  to  apply  to  the  case  of  van¬ 


ishingly  small 


4.  2  Forward  Flight  Theoi 


As  in  the  hovering  case  the  governing  equations 


i  =■  -  5 


>5  ...  i*  ,  j 


(60) 


(61) 


where 


now  we  put 


\JL. 

»  v'  )X 


(62) 


Once  again  these  equations  are  linearised,  but 


i,  -  <af  (l 
v;  » 


(63) 


(64) 


where  is  the  tip  speed  ratio  and  is  the 

mean  inflow  velocity  coefficient.  Putting  equations  (60)  - 
(65)  in  non-dimensional  form  by  means  of  the  transformations 
(38)  ,  equation  (60)  becomes 


31 


(66) 


and  equation  (61)  becomes 

$  +-0tA*'n^Tj  +  ^  7^  -  ° 

To  solve  (66)  we  expand  in  ascending  powers 
of  JJL  i.e.  ,  we  put 

<f>-  (67) 

This  is  a  common  procedure  in  the  treatment  of  differential 

/16) 

equations  with  periodic  coefficients'1  ,  and  it  has  many 
applications  to  the  helicopter  in  forward  flight  (see  Jones 
and  Shutler^^). 


Substituting  (67)  in  (66)  and  writing  the  whole 
as  an  expansion  in  powers  of  l/i  we  obtain 

4-sin  ty)U,  -  o 

If  (68)  is  to  be  valid  for  all  values  of  then  the 

coefficient  of  each  power  of  jL  must  be  separately  zero. 
Thus  (68)  becomes  the  following  set  of  successive  partial 
differential  equations 


(68) 


plus  etc.,  for  higher  powers  of  A.  than  the  first.  For 
the  rest  of  this  analysis  only  the  first  power  of  will 


(69) 

(70) 
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be  retained.  This  should  be  sufficient  to  give  the  rate 
of  change  with jL  as Ji.  tends  to  zero,  and  the  whole  approH- 
imation  does  not  warrant  any  further  complication. 

It  will  be  seen  that  the  equation  for  ^is  iden¬ 
tical  with  the  equation  obtained  from  the  vorticity  distri¬ 
bution  in  the  hovering  case.  The  difference  is  that  in 
forward  flight  many  frequencies  will  be  present  in  the  in¬ 
duced  velocity  and  blade  motion  and  these  will,  of  course, 
appear  in  the  vorticity  distribution. 


Hence  we  put 


and  at  the  same  time 

Kcot't 

i  - 1  S'  V  *•» .*f + 1  £  2v>ec#SY>  f 

The  substitution  (71)  makes  ^independent  of  ^  »  which  is 
as  it  must  he  since  the  inter-blade  phase  relationship  is 
such  as  to  make  for  each  frequency. 


(71) 


(72) 


becomes 


Substituting  (71)  in  (70)  the  equation  for 
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Thus  to  a  first  order  in  A 

4  *  (H)  +MT«  f  sin  n(f-j)  ] 

4_  C*  £  £*>«•*>  (f-  j)  -yOl  cos^cm  n  -T>J 

(75)  is  now  substituted  in  the  partial  differen¬ 
tial  r  ion  (65)  to  give  the  non-dimensional  stream  func- 


(75) 


tion 


solution  for  F  is 

Z. 


Since  (75)  is  independent  of 


X 


the  complete 


OSS  f  ■‘Si*  V1 


(Mi 

*  _  \IVI  J  . 


(76) 


+-X  C0  e,1t«'((,l)+S^,e 

wher^  the  complementary  function  has  been  chosen  to  ensure 
finite  velocities  at  infinity. 

Since  there  is  no  vorticity  above  the  disc 

f  *  "(H) + r ’ V*7* " "  (H) 

'  V  * 

The  equation  of  continuity  at  the  actuator  disc  is 

>5  If 

which  when  evaluated  at  the  reference  blade  becomes 


(77) 


(78) 
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Y  Is  £>,®,  sin  n  fj- 

r»  _  _  (79) 

$» c**"  4-^  V*  ^s*rc'»s"'  4 

^  Similarly  the  equation  (32)  which  expresses  the 

condition  that  changes  in  circulation  around  a  blade  are 
convccted  away  normal  to  the  disc  becomes 

T^fen+O  <>H  4  +%**(S*+$d  f  (80) 

Finally  the  equation  (37)  connecting  the  lifting 
line  and  actuator  disc  theories  for  the  non-steady  part  of 
the  motion  becomes 


("V  CM*/  +£*  *'•*  ^ 

■n,  * 

•# 


(81) 


cas  h 
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In  each  of  equations  (79),  (80),  (81)  the  various 
frequencies  are  collected  together,  and  the  coefficients  of 

are  equated  separately  to  zero  to  gi' 
equations  for  }  § 


-ve  six 


ns 


in  terms  of  )  etc.  Eliminating  the  arbitrary 

constants  the  following  equations  are  obtained 

v-  l  V  J  T  ^  J 


_  z. 

V 


(82; 
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^ss'+  ^jjK-04+^3 

In  arriving  at  (82)  ,  (83)  terms  have  been  neglected 

on  the  grounds  that  A  and  ^  are  of  the  same  order,  and 
terms  in  r  have  already  been  ignored  as  too  small. 


terms  in 


have  already  been  ignored  as  too  small. 


If  we 


put  A*?  in  (82)  ,  (83)  the  equa- 

iquations  (53)  ,  (54)  for  the  vorticity 


since 


14-vicjy^ 


depends  upon  the  quantities  i)  i.e.,  upon  the 
circulation  at  adjacent  frequency  ratios.  This  particular 
harmonic  coupling  has  it  origins  in  the  varying  convection 
speed  of  the  vorticity.  When  the  lift  is  calculated  there 
is  an  additional  inter-harmonic  coupling  since 


l_=ry0w) 


(83) 


tions  reduce  to  equations  (53)  ,  (54)  for  the  vorticity 
components  in  the  hover.  To  0  (jty  therefore  the  circula¬ 
tion  at  a  particular  frequency  ratio  which  derives  from 


(84) 


(85) 


Therefore,  whereas  in  hovering  it  is  possible  to  relate  the 
lift  at  a  given  frequency  ratio  directly  to  the  blade  motion 
at  that  frequency  ratio,  in  forward  flight  this  cannot  be 
done.  Instead  it  is  necessary  to  carry  the  calculation  one 
stage  further  and  solve  the  dynamical  equations  of  motion 
for  the  blades,  taking  into  account  a  whole  spectrum  of 
modes  and  harmonics. 
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V 


DISCUSSION 


A  comparison  of  Equations  (82)  and  (83)  with  equa¬ 
tions  (53) ,  (54)  shows  that  forward  flight  affects  the  shed 
vorticity  -  and  hence  the  circulation  around  the  blades  - 
directly  through  the  induced  velocity  variations,  but  it 
also  introduces  a  cross-coupling.  Thus,  for  example,  the 
third  harmonic  components  of  circulation  around  a  blade  are 
determined  by  the  third  harmonic  components  of  the  relative 
induced  velocity,  together  with  some  contribution  from  the 
second  and  fourth  harmonic  circulation  components.  These 
inter-harmonic  couplings  are  indeed  proportional  to  JL.  but 
a  closer  inspection  of  Equations  (82)  and  (83)  shows  their 
effects  to  depend  upon  the  much  larger  quantity •  It 
seems  unlikely  therefore  that  the  effects  of  inter-harmonic 
coupling  can  be  ignored  although  it  may  well  be  that  this 
model  of  the  flow  inflates  their  influence.  On  the  other 
hand,  the  inter-harmonic  coupling  in  the  lift  which  arises 
from  the  interaction  between  the  fluctuations  in  relative 
wind  speed  and  the  variations  in  the  circulation  -  sec 
equation  (85)  -  is  only  proportional  to 


A* 


The  numerical  significance  of  the  couplings  is 
not  easy  to  see  because  the  individual  circulation  components 
should  really  be  obtained  by  solving  a  large  number  of  sim¬ 
ultaneous  equations.  But  some  special  cases  and  features 
can  be  examined  and  these  are  most  illuminating. 


In  lifting  line  theory  is  employed,  the  blade 
airloads  at  a  particular  frequency  ratio  are  proportional 
to  the  right-hand  sides  of  (82),  (83).  After  solving  the 
equations  the  constant  of  proportionality  is  found  to  be 
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the  effective  lift  curve  slope  or  real  part  of  the  lift  de¬ 
ficiency  function.  For  rotor  bl.ides  this  is  usually  regarded 
is  being  a  function  only  of  A  and.  ^  .  But  if  the  inter¬ 

harmonic  coupling  is  large,  equations  (o2)  and  (33)  show  that 
the  constant  of  proportionality  is  no  longer  simple.  It  must 
depend  upon  the  circulation  content  at  other  harmonics,  which 
in  turn  depend  upon  the  blade  motion  and  the  induced  velo¬ 
city  distribution  at  these  harmonics.  Our  first  conclusion 


then  is  that,  within  the  limits  of  velocity  of  the  present 
theory,  it  is  not  legitimate  to  isolate  the  lift  curve 


slope  from  the  blade  motion  and  induced  velocity.  More 


precisely 


harmonics  in 


can  only  be  defined  if  the  proportion  of 
is  specified. 


The  second  conclusion  we  can  draw  concerns  the 
damping  cf  blade  motion.  Since  Equations  (82)  and  (83)  are 
linear,  we  can  consider  the  effect  of  a  small  disturbance 
in  the  blade  motion  separately.  If  this  is  confined  to  a 
single  harmonic,  say  the  so  that  )  etC *  ’ 

are  all  zero  then,  because  of  the  cross-coupling,  the  cir¬ 
culation  at  the  representative  section  will  contain  all 
harmonics.  The  equations  for  these  additional  harmonics 
will  be  of  the  type 

jb  0}  hiF 


The  proportion  of  will  depend  upoi 

A  and  the  magnitude  of  thej^  ^  t^lc 

theory  is  to  be  valid  the  etc.,  must  be  insig¬ 

nificant.  Now  suppose  the  disturbed  blade  motion  to  take 


will  depend  upon 


but  if  the 


(86) 


(87) 
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place  at  a  frequency  ratio  i.e.,  the  natural  frequency 

of  the  blade  is  close  to  cycles  per  revolution.  This 
means  that  the  generalised  inertia  force  at  this  frequency 
almost  balances  the  generalised  elastic  force,  hence  the 
generalised  aerodynamic  force  must  be  small.  Therefore  in 
the  present  model  must  6e  small  and,  from  (86),  (87) 

will  be  even  smaller.  Very  roughly 


^1 


n* 


(88) 


(89) 


Substituting  (88) ,  (89)  and  the  corresponding 
expressions  for  .  v  in  to  (82)  ,  (83)  shows  that  the 

'correction*  to  the  interharmonic  coupling  will  be 

of  order  (m*0z  and  therefore  negligible.  Equations  (82) 
and  (83)  then  become  identical  with  the  hovering  equations 
(53) ,  (54)  which  lead  directly  to  the  result  C?=ly/k+Tr  for  the 
lift  deficiency  function.  We  conclude,  therefore,  that 
the  damping  of  natural  oscillations  at  a  particular  fre¬ 
quency  is  unaffected  by  UL  for  small 


This  result  may  seem  rather  surprising,  but  there 

is  some  experimental  evidence  to  support  it,  and  it  can  be 

justified  theoretically  in  another  way.  The  experimental 

(18) 

evidence  is  due  to  Ham,  Moser  and  Zvarav  ' .  They  measured 
the  flapping  response  of  the  blades  of  a  model  helicopter 
rotor  to  a  vertical  excitation  of  the  hub.  In  the  hovering 
case,  for  low  blade  pitch  angles,  this  response  showed,  in 
accordance  with  'wake*  theory,  a  sharp  peak  at  a  frequency 
of  two  cycles  per  revolution  (the  rotor  was  of  the  two-bladed, 
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hinged  type) .  In  forward  flight  this  peak  persisted  even 

Theoretical  justification*  can  be  got  by  arguing 
that  the  flow  model  and  the  real  rotor  system  are  reversible, 
i.e.,  it  does  not  matter  whether  U  is  positive  or  negative. 
Hence  a  plot  of  the  real  part  or  the  Theodorsen  function 

.(14), 


inst  U  ,  such  as  that  given  by  Miller v'L^/must  be  symmetric 
ut  the  axis  ,  i.e.  (UiM  =  0  and  the  change  with 


This  result  is  sketched 


aga 
about 

A  can  at  most  be  of  order  A  . 

m  Figure  3,  superimposed  upon  the  results  of  Miller  s  cal¬ 
culation.  Of  course  the  sketch  is  conjecture  -  except 

for  the  slope  at  0  -  but  the  horizontal  part  has  been 

projected  to since  the  limit  of  this  theory  is 


Away  from  a  natural  frequency  though  this  result 
may  not  hold;  there  is  in  fact  a  difficulty  in  the  concept 
of  damping  in  such  cases.  Usually  damping  is  taken  to  mean 
the  coefficient  of  a  velocity  -  in  this  case  at  a  particular 
frequency  ratio  or  harmonic  Vr\ .  This  quantity  may  perhaps 
be  unchanged  by  forward  motion  but  an  isolated  harmonic  can¬ 
not  exist.  The  inter-harmonic  coupling  will  ensure  that 
blade  motion  and  or  vortex  shedding  will  occur  at  other  fre¬ 
quencies.  It  will  be  necessary  to  do  work  on  the  system 
to  maintain  these  motions  and  manufacture  the  vorticity. 
Hence  a  more  appropriate  measure  of  the  damping  would  prob¬ 
ably  be  the  energy  dissipated  in  maintaining  an  oscillation 
and  this  will  certainly  increase  with  • 

When  we  come  to  consider  the  more  important  prob¬ 
lem  of  forced  oscillation  due  to  the  induced  velocity  var¬ 
iations,  we  have  to  recognise  that  resonance  is  difficult 

*1  am  indebted  to  Dr.  Statle  of  the  Cornell  Aeronautical 
Laboratory  for  this  explanation. 
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to  avoid  in  practice.  Now  the  generalised  aerodynamic  force 
is  zero  at  resonance ,  and  if  there  is  no  inter -harmonic 
coupling  tne  blade  incidence  must  be  effectively  zero  at 
the  resonant  frequency.  This  means  that  the  blade  bends 
in  sympathy  with  the  induced  velocity  variations  so  that 
there  are  no  net  incidence  changes  at  that  frequency.  But 
if  the  inter-harmonic  coupling  is  strong,  then  without  solv¬ 
ing  the  equations  wc  cannot  say  what  the  blade  motion  will 
be.  The  circulation  component  at  the  resonant  frequency 
will  be  small,  probably  much  less  than  at  the  adjacent  fre¬ 
quencies  (since  they  will  not  be  influenced  by  resonance). 
Therefore  it  is  likely  that  the  inter -harmonic  coupling 

terms  will  determine  the  blade  motion  at  resonance.  This 

(14) 

possibility  was  forseen  by  Miller'  . 
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VI 


CONCLUSIONS:  FURTHER  DEVELOPMENTS 


A  comparatively  crude  model  of  the  flow  through 
a  rotor  has  been  developed  for  estimating  the  influence  of 
the  shed  wake  at  low  tip-speed  ratios.  It  gives  results 
which  are  in  agreement  with  more  refined  theories  for  the 
hovering  case,  and  there  are  reasonable  grounds  for  expect¬ 
ing  the  model  to  be  valid  at  lease  as jJ^Q»  Therefore,  al¬ 
though  this  particular  application  leaves  something  to  the 
imagination,  it  is  not  felt  that  any  further  apology  or 
just:'  fication  for  its  use  is  necessary.  A  well-established 
result  from  hovering  theory  has  been  reproduced,  and  there 
is  adequate  physical  explanation  for  many  of  the  theoretical 
features  which  are  brought  out.  For  instance  the  parameter 
appears  quite  clearly  as  a  measure  of  how  much  the  wake  is 
'filled'  by  shed  vorticity,  the  fact  that  can 

be  deduced  from  quite  separate  considerations  of  symmetry, 
and  even  the  dependence  of  the  inter-harmonic  coupling  on 
the  product  of  tip  speed  ratio  and  harmonic  number  can  be 
explained. 


The  extension  to  forward  flight  proceeds  by  an 
expansion  in  powers  of  but  analysis  shows  that  the  ex¬ 
pansion  really  should  be  in  powers  of  yOn. ;  tvl  here  is  the 
(integral)  order  of  vibration  or  frequency  ratio.  Because 
of  this,  it  is  found  that  the  inter-harmonic  couplings, 
i.e.,  the  generation  of  aerodynamic  forces  at  one  frequency 
by  circulation  changes  at  another,  are  surprisingly  large. 
Their  influence  is  strong  enough  to  suggest  that  the  simple 
idea  of  a  lift  curve  slope  independent  of  the  details  of 
the  blade  motion  and  induced  velocity  variati  should 

be  abandoned.  In  particular,  the  interharmo  plings 


must  be  taken  into  account  if  one  or  more  blade  natural  fre¬ 
quencies  lie  close  to  an  integral  number  of  cycles  per  revo¬ 
lution. 


For  oscillations  at  a  single  frequency  at  or  close 
to  a  natural  frequency,  the  effect  of  interharmonic  couplings 
can  be  ignored.  It  is  then  found  that  the  damping  of  the  os¬ 
cillations,  defined  as  the  aerodynamic  coefficient  of  a  blade 
bending  velocity  at  that  frequency,  is  unaffected  by  . 

This  may  account  for  the  apparent  persistence  of  shed  wake 
effects  as  the  tip  speed  ration  increases. 

But  care  is  necessary  in  interpreting  damping  in 
this  case,  since  the  oscillation  at  a  single  frequency  is 
accompanied  by  vortex  shedding  at  the  coupled  frequencies. 
Work  will  undoubtedly  have  to  be  done  on  the  flow  to  pro¬ 
duce  the  shed  vorticity  and  this  will  appear  as  an  energy 
dissipation  somewhere  else.  In  wind  tunnel  experiments, 
where  model  rotors  are  forced  to  oscillate  b  some  external 
agency,  there  will  be  three  possible  sources  of  energy  viz. 
the  tunnel  fan  motor,  the  rotor  drive  and  the  exciter  power 
supply.  Hence  it  may  not  be  strictly  valid  to  attribute 
large  oscillations  to  small  aerodynamic  damping  unless 
measurements  of  the  total  power  consumption  are  made. 

It  remains  now  to  exploit  this  concept  further, 
either  by  straightforward  application  of  the  actuator  disc 
to  other  problems  or,  perhaps  better,  by  an  extension  of 
the  use  of  momentum  and  energy  methods.  These  cannot  be 
expected  to  yield  detailed  answers  but  they  do  highlight 
those  features  which  a  more  elaborate  calculation  should 
include.  In  the  spectrum  of  techniques  vh  ich  are  generally 
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available:  they  lie  somewhere  between  dimensional  analysis 
and  numerical  solution,  and  as  such  they  should  provide  a 
gcod  guide  to  future  exocrirnental  work. 
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Fig.  1  -  Nomenclature  for  flow  model 


Fig.  2  -  Explanation  for  forward  flight  model 


Fig.  3  -  Variation  of  F  with  tip  speed  ratio 
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